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摘 要 : 客观 地 评价 水 文 情势 变化 是 进行 河流 生态 修复 和 流域 水 资源 综合 管理 的 基础 。 改 进 RVA 法 将 径流 典型 


年 份 的 变化 纳入 整体 改变 度 的 计算 过 程 中 ,相对 更 全 画 


ij 地 评价 径流 水 文 情势 的 变化 。 近 几 十 年 来 黑河 流域 气候 条 


件 的 发 生 显著 变化 ,黑河 上 游 河流 水 文 情势 的 变化 情况 仍 需 进一步 研究 。 以 黑河 上 游 为 研究 区 域 ,应 用 M-K 检 验 


识别 径流 的 突变 年 份 ,采用 改进 RVA 法 计算 分 析 黑 河上 游 河 流水 文 情势 变化 。 结 果 表 明 ,黑河 上 游 径流 的 突变 年 


份 为 2002 年 ,2002 年 前 后 分 别 以 平水 年 和 丰 水 年 为 主 ; 采 
法 计算 的 整体 改变 度 为 74.62% ,说 明 RVA 法 低估 了 黑河 上 游 水 文 情势 整体 改变 度 ; 黑 河上 游 径流 水 文 情势 变化 主 
要 包括 : 相 比 于 1960 一 2002 年 ,2003 一 2015 年 后 各 月 平均 流量 均 呈 不 同 程度 增加 ,但 最 小 1 日 3 日 .7 日 和 最 大 7 日 
流量 减 小 ,其 余 极 端 流量 指标 呈 不 同 程度 增加 ,高低 流 量 发 生 次 数 及 历时 变化 微弱 ,流量 平均 增加 率 , 减 小 率 以 及 


i 


J RV A 法 计算 的 整体 改变 度 为 54.68% ,而 改进 的 RAV 方 


逆转 次 数 均 增加 ;采用 改进 的 RVA 法 计算 的 黑河 上 游 水 文 情势 的 整体 改变 度 为 74.62% ,属于 高 度 改变 。 气 候 变 化 


是 导致 黑河 上 游 水 文 情势 高 度 改变 的 主要 原因 。 


关键 词 : 河流 水 文 情势 变化 ; RVA 法 ; 典型 年 份 ; 黑河 上 游 


水 文 情 势 决定 并 影响 着 河流 生态 系统 的 物质 
循环 .能 量 过 程 物理 栖息 地 状况 。 随 着 气候 变化 
效应 凸显 ,人 类 活动 加 剧 ,河流 的 水 文 情势 发 生 了 
显著 变化 ,这 些 变 化 直接 影响 河流 生态 系统 的 生物 
群落 组 成 及 结构 ""。 客 观 地 评价 水 文 情势 变化 是 进 
行 河流 生态 修复 和 流域 水 资源 综合 管理 的 基础 。 
评价 河流 水 文 情势 变化 的 常用 方法 为 水 文 变化 指 
brid (Index of Hydrological Alteration , IHA ) 以 及 在 此 
基础 上 发 展 的 变动 范围 法 (Range of Variability Ap- 
proach, RVA 法 )。 目 前 ,RVA 法 已 经 广泛 应 用 到 
内 外 河流 径流 情势 变化 评价 中 。Galat 等 ”对比 了 
1929—1948 年 与 1967 一 1996 年 两 个 时 期 密苏里 河 
11 个 观测 站 点 径流 水 文 情 势 的 变化 。Shiau 等 “ 评 


的 影响 。 李 兴 拼 等 研究 了 枫 树 坝 水 库 的 运行 对 东 
江上 游 径流 情势 的 影响 ,结果 表明 东江 上 游 的 径流 
AAA EME, RES DPSS” yA 
RVA 法 对 河流 水 文 情 势 的 变化 进行 了 分 析 。 相 关 
研究 表明 ,RVA 在 计算 整体 改变 度 时 将 变化 前 后 两 
个 时 间 段 各 视 为 一 个 整体 ,忽略 了 径流 典型 年 份 的 
变化 ,会 导致 对 水 文 情势 变化 判断 的 不 准确 ,进而 
对 水 资源 管理 产生 不 利 影响 ”。 本 文采 用 改进 的 
RVA 法 ,将 径流 典型 年 份 的 变化 纳入 整体 改变 度 的 
计算 过 程 中 ,可 相对 更 全 面 地 评价 径流 水 文 情势 的 
变化 。 

黑河 是 我 国 第 二 大 内 陆 河流 ,中 游 张 掖 地 区 是 
古 丝 绸 之 路 . 今 欧 亚 大 陆桥 要 地 ,下游 额济纳 地 区 


价 了 分 水 卉 的 修建 对 Chou-Shui Creak 河 的 影响 。 
Irwin 等 外 评价 了 1924—1980 年 和 1984 一 1996 年 
Tallapoosa 河流 量 特征 的 变化 。 郭 文献 等 中 采用 
RVA 法 分 析 了 丹江口 水 库 对 汉江 中 下 游 径流 情势 
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是 我 国 重要 的 国防 科研 基地 。 流 域内 生态 环境 建 
设 与 保护 还 关系 到 西北 .华北 地 区 的 环境 质量 , 具 
有 重要 的 战略 地 位 。 黑 河流 域 水 资源 时 空 分 布 严 
重 不 均 ,流域 水 资源 供需 矛盾 突出 ,是 限制 该 区 
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经 济 社会 发 展 的 主要 障碍 ,也 是 造成 其 生态 环境 极 
度 脆弱 的 主要 原因 。 黑 河流 域 70% 左 右 的 径流 产 
生 于 营 落 峡 水 文 站 以 上 的 上 游 地 区 , 且 中 下 游 90% 
以 上 的 水 资源 都 依赖 上 游 补给 "”。 黑 河上 游 的 水 
文 情势 变化 对 整个 流域 水 资源 具有 重要 影响 。 

国内 外 学 者 对 黑河 流域 水 文 特征 及 生态 环境 
变化 开展 了 一 些 研究 。 王 根 绪 等 "对 比分 析 了 20 
世纪 40 年 代 以 来 不 同时 期 黑河 流域 水 文 及 生态 环 
境 的 变化 特征 ,并 对 不 同时 期 流域 水 土 资源 开发 规 
模 及 出 山 径流 的 变化 进行 了 分 析 ,结果 表明 ,黑河 
流域 水 系 发 生根 本 性 变化 。 李 弘毅 等 ” 综 合 研究 
了 黑河 上 游 祁 连 山 冰 沟 流域 积 雪 水 文 情势 ,并 分 析 
了 冻 土 水 热 变 化 过 程 和 融雪 径流 变化 特征 。 郭 巧 
玲 等 ”研究 了 黑河 流域 径流 年 内 分 配 变化 。 党 素 
珍 等 "研究 了 黑河 上 游 的 基 流 及 其 变化 特征 ,表明 
黑河 上 游 及 流量 年 内 变化 先 增 后 减 ,年 际 变化 呈 增 
加 趋势 ,而 基 流 指数 呈 减 少 趋势 。 邹 悦 等 ”研究 了 
黑河 流域 营 落 峡 站 河道 内 生态 需 水 ,指出 1997 年 前 
后 河道 内 生态 需 水 和 满足 生态 需 水 出 现 频率 变化 
显著 。 程 国 栋 等 ”对 黑河 流域 生态 -水 文 过 程 集成 
开展 了 研究 ,初步 揭示 了 流域 生态 水 文 过 程 耦合 机 
理 ,形成 了 黑河 关键 地 区 水 循环 系统 的 基本 认识 。 
Sang 等 "分析 了 变化 环境 对 黑河 上 中 游 径流 的 影 
响 , 指 出 上 游 降水 增加 和 温度 的 上 升 导致 径流 的 增 
加 。 王 宇 涵 等 "研究 了 冰冻 圈 水 文 过 程 对 黑河 上 
游 径流 的 影响 ,指出 黑河 上 游 在 近 54 a 内 呈现 上 升 
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趋势 ,其 中 以 基 流 上 升 为 主 。 李 秋 菊 等 ”分析 了 黑 
河上 游 径流 变化 及 其 原因 ,并 量化 了 降水 ` 冰 川 融 
雪 以 及 蒸 散 发 对 径流 增加 的 贡献 。 随 着 近 儿 十 年 
来 黑河 气候 条 件 的 显著 变化 ,上 游 河流 水 文 情势 的 
变化 情况 仍 需 进 一 步 研 究 。 本 文 以 黑河 上 游 为 研 
究 区 域 ,收集 营 落 峡 水 文 站 1960 一 2015 年 逐日 观测 
流量 ,应 用 Mann-Kendall 检验 (M-K 检 验 ) 识 别 径流 
的 突变 年 份 ,采用 改进 RVA 法 计算 了 水 文 情势 的 整 
体 改变 度 ,并 分 析 了 其 水 文 情 势 变化 的 主要 原因 ， 
为 流域 河流 生态 修复 和 水 资源 综合 管理 提供 参考 。 


1 研究 区 概况 及 数据 来 源 


1.1 研究 区 概况 

本 文选 取 黑 河 干 流 上 游 为 研究 区 ,其 位 置 及 
DEM 如 图 1 所 示 。 黑 河 干流 以 营 落 峡 和 正义 峡 为 界 ， 
分 为 上 游 .中游 和 下 游 。 营 落 峡 水 文 站 以 上 为 上 游 ， 
干流 上 游 面积 约 10000 km ,海拔 高 度 约 在 1600 ~ 
5100 m。 该 地 区 气候 寒冷 ,年 平均 温度 为 -5~4°%C""。 
土地 利用 类 型 以 森林 草地 岩石 和 裸 地 为 主 。 
1.2 数据 来 源 

收集 营 落 峡 水 文 站 1960 一 2015 年 逐日 观测 流 
量 数据 以 及 研究 区 周边 4 个 气象 站 的 逐日 降水 数 
据 , 研 究 区 水 文 站 与 气象 站 的 位 置 如 图 1 所 示 。 数 
据 源 于 “国家 青藏 高 原 科学 数据 中 心 ”(http://westde. 


westgis.ac.cn/) . 
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图 1 研究 区 位 置 


Fig. 1 Location and DEM ofthe study area 
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2 研究 方法 


黑河 上 游 水 文 情势 变化 分 析 主 要 分 为 以 下 几 
个 步骤 :首先 ,采用 M-K 检 验 识别 径流 的 突变 年 份 ， 
将 观测 数据 划分 为 突变 前 后 两 个 水 文 序列 ;然后 ， 
采用 改进 的 RVA 方 法 计算 河流 水 文 情势 的 整体 改 
变 度 ,具体 步骤 包括 :选取 IHA 法 分 析 河 流水 文 情 
势 的 变化 情况 ,该 方法 包含 反映 数量 大 小 发生 时 
间 持续 时 间 、 频 率 以 及 变化 率 等 特征 的 5 组 共计 33 
个 水 文 指标 ;在 IHA 法 基础 上 ,采用 RVA 法 计算 33 
个 水 文 指 标的 改变 度 ; 采 用 欧 拉 距离 计算 径流 突变 
年 份 前 后 典型 年 份 ( 丰 平 . 枯 水 年 ) 的 变化 ,结合 
IHA 法 计算 的 水 文 指 标的 改变 度 和 欧 拉 距离 ,计算 
河流 水 文 情势 的 整体 改变 度 。 
2.1 M-K 检 验 

M-K 检 验 最 初 由 Mann FU Kendall fE, ZN 
法 在 时 间 序 列 趋 势 分 析 及 突变 点 识别 时 ,不 要 求 样 
本 遵循 一 定 的 分 布 ,也 不 受 异 常 值 的 干扰 且 计 算 简 
便 , 因 此 ,受到 世界 气象 组 织 推荐 ,普遍 应 用 于 水 
文 .气象 等 非 正 态 分 布 的 数据 分 析 失 2” 。 本 文选 用 
该 方法 进行 突变 点 的 识别 。 
2.2 改进 的 RVA 法 

改进 RVA 方 法 首先 由 Yin 等 提出 ,在 常用 的 
河流 水 文 情 势 变化 评价 方法 RVA2 2 的 基础 上 考虑 
了 径流 典型 年 份 的 变化 ,使 计算 结果 相对 更 合理 。 
该 方法 通过 欧 拉 距离 计算 典型 年 份 的 变化 ,再 将 
计算 结果 代入 水 文 情势 改变 度 的 计算 中 ,从 而 将 径 
流 典 型 年 份 的 变化 纳入 水 文 情势 变化 评价 当中 。 
改进 RVA 方 法 具体 步骤 及 计算 如 下 : 

丰 平 .村 水 年 分 别 赋 值 2、1、0。 径 流 突变 前 后 
典型 年 份 的 分 布 情况 为 (41,4;,…,A,) 和 (B,B,,…， 
B,) ,A;: 和 B=0、1、2, 欧 拉 距 离 (ED) 的 公式 ”为 : 


如 果 a=b， 那 么 ，ED= YB) (1) 


WRa>b, WBA, ED=min Bear , 


k=0, 1, 2:.n—k 
(2) 


如 果 a<b， IA, ED=min X (4,-B Y ， 


k=0, 1, 2:--n-k 
(3) 


RP: HAM BEEKE .平水 年 或 枯 水 年 时 ， 
|A,-B,|=0 ; MAMB P—TRAPPAE, AT AE 
水 年 或 者 枯 水 年 时 , |A,-B,]=1 5; 4A ANB TSI AT 
水 年 和 丰 水 年 或 者 丰 水 年 和 枯 水 年 时 ， 
|4,-B,|=2. 

将 计算 出 的 欧 拉 距离 进行 标准 化 ,公式 如 下 : 

S= EB (4) 
Tx min (va, Vb) 
式 中 :5 表示 标准 化 后 的 欧 氏 距离 ,用 于 衡量 径流 的 
丰 . 平 、 枯 水 年 的 变化 ,变化 范围 为 L0,1j ,数值 越 大 
表示 变化 越 大 ;ED 表示 欧 拉 距离 ;7 表示 典型 年 份 
的 最 大 距离 , 当 有 3 种 典型 年 份 ( 丰 、 平 .村 水 年 ) 时 ， 
T=2, 

Grantz 等 2 定义 丰 水 年 为 流量 大 于 或 等 于 总 记 
录 第 75 百 分 位 流量 的 年 份 , 枯 水 年 为 流量 小 于 总 记 
录 第 25 百 分 位 流量 的 年 份 ,流量 介 于 两 者 之 间 的 年 
份 为 平水 年 。 文 草 采 用 该 方法 ,将 1960 一 2015 年 共 
计 56 a 划 分 为 丰 水 年 .平水 年 和 村 水 年 。 

考虑 典型 年 份 变 化 的 整体 改变 度 的 计算 公式 
A”: 


OD=1-(1-D,)d-S) (5) 
式 中 ;0D 表示 整体 改变 度 ,0D 值 介 于 0~33% 时 为 低 
度 改 变 , 介 于 33%~67% 时 为 中 度 改变 , 介 于 67%~ 
100% 为 高 度 改变 ;Do 表示 RVA 法 计算 的 整体 改变 
度 ;S 表 示 标 准 化 后 的 欧 氏 距离 。 
Do 具体 计算 公式 为 : 


4 
Di= 丁 六 D， (6) 
m=1 


式 中 :4 表示 IHA 水 文 指标 总 数 ,IHA 水 文 指标 主要 
包括 月 平均 流量 .年 极端 流量 .年 极端 值 发 生 时 间 、 
高 低 流 量 发 生 频 率 及 持续 时 间 和 流量 变化 改变 率 
及 频率 (次 )j5 大 类 , 共 33 个 指标 2 5 D, BEAN SB m 
IHA 指 标的 改变 度 。 
每 个 IHA 水文 指 标 改 变 度 的 计算 公式 "中 如 下 : 
N,-N, 
N 


D,,= x 100% (7) 


N,=r,,XN, (8) 
式 中 :D, 表 示 第 mn 个 水 文 指 标的 改变 度 ,D, 值 介 于 
0~33% 时 为 低 度 改变 , 介 于 33%~67% 时 为 中 度 改 
变 , 介 于 67%~100% 为 高 度 改变 ; 入 ,表示 第 m 个 指标 
在 变化 后 实际 落 入 底 化 范围 内 的 年 数 ,NN. 表 示 第 m 
个 指标 变化 后 预期 落 入 变动 范围 内 的 年 数 ;mm 表 示 
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变化 前 第 mm 个 指标 落 在 变化 范围 内 的 比例 , 即 变 化 
前 第 m 个 指标 落 入 变化 范围 内 年 数 与 变化 前 流量 
记录 总 年 数 的 比值 ;入 , 表 示 变 化 后 流量 记录 的 总 
年 数 。 


3 结果 与 分 析 


3.1 黑河 上 游 水 文 指标 变化 

3.1.1 黑河 上 游 年 均 流 量 突变 点 的 识别 FH M-K 
检验 方法 对 营 落 峡 水 文 站 1960 一 2015 年 均 流 量 数 
据 进 行 突变 检验 ,UF 和 UB 曲线 如 图 2 所 示 。 结 果 
表明 ,UF 与 V8 仅 有 一 个 交点 ,该 交点 的 横 坐 标 为 
2002 年 ,表明 营 落 峡 站 流量 在 2002 年 发 生 了 突变 。 
因此 ,将 1960—2015 年 分 为 1960—2002 年 和 
2003 一 2015 年 两 个 时 期 。1960 一 2015 年 均 流量 如 
图 3 所 示 , 1960—2002 年 均 流 量 为 48.6 ms ,而 


统计 量 


1970 1980 


1990 
年 份 


2003—2015 年 均 流 量 为 60.5 ms ,观测 流量 增幅 
H 24.49%. 
3.1.2 突变 前 后 的 水 文 指标 的 变化 “基于 M-K 方 法 
识别 的 突变 点 将 观测 径流 数据 划分 为 两 个 水 文 序 
列 ,并 对 两 个 水 文 序列 采用 RVA 进行 分 析 ( 表 1) ,各 
旧 标 改变 程度 如 图 4 所 示 。 

(1) 各 月 平均 流量 。 由 表 1 可 以 看 出 ,突变 后 
营 落 峡 水 文 站 各 月 的 平均 流量 均 增加 。 其 中 ,1 月 、 
2 月 和 5 月 流量 的 变化 率 较 小 ,接近 10%;8 一 11 月 


流量 的 变化 率 最 大 ,在 50% 左 右 , 表 明 8 一 11 月 平均 
流量 的 增加 最 为 明显 。 各 月 平均 流量 改变 度 中 ,5 
月 .6 月 最 小 , 均 为 15.05%,7 月 份 为 28.09%, 属 于 低 
度 改 变 ;1 月 .2 月 、12 月 的 改变 度 介 于 33% ~ 67%, 
属于 中 度 改 变 ;其 余 月 份 的 改变 度 均 大 于 67%, 属 
于 高 度 改 变 。 


2000 2010 2020 


图 2 黑河 上 游 营 落 峡 水 文 站 流量 的 M-K 检 验 
Fig. 2 M-K test of the observed runoff at Yingluoxia hydrologic station 


100 


80 


60 48.62 m's! 


日 均 流量 /(m3s-1!) 


40 


1960 


1980 


1970 


图 3 黑河 上 游 营 落 峡 水 文 站 1960 一 2015 年 日 均 流量 
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Fig. 3 The mean annual flow at Yingluoxia hydrologic station from 1960 to 2015 
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表 1 萤 落 峡 水 文 站 观测 流量 IHA 水 文 指标 统计 表 


Tab. 1 Statistics of IHA index at Yingluoxia hydrologic station 
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THA 指标 类 别 水 文 指标 变化 前 变化 后 变化 率 /% 改变 度 /% 程度 
月 平均 流量 1 月 12.9 14 8.53 42.47 中 度 
2 月 13.69 14.65 7.01 42.47 中 度 
3 月 15.59 22.6 44.96 71.24 高 度 
4 月 22.88 32.15 40.52 85.62 高 度 
5 月 37.8 41.2 9.00 15.05 低 度 
6 月 61.45 74.05 20.51 15.05 低 度 
7 月 103.9 131 26.07 28.09 低 度 
8 月 87 129 48.28 71.24 高 度 
9 月 69.35 102.8 48.16 71.24 高 度 
10H 35.27 52.6 49.14 100 高 度 
1H 22.5 34 51.11 100 高 度 
12H 14.9 16.9 13.42 58.65 中 度 
年 极端 流量 最 小 1 日 8.08 0 -100.00 85.62 高 度 
最 小 3 日 8.757 0 -100.00 85.62 高 度 
最 小 7 日 9.473 4.187 -55.80 85.62 高 度 
最 小 30 日 11.74 11.89 1.26 28.09 低 度 
最 小 90 日 14 16.81 20.09 56.86 中 度 
最 大 1 日 351 404 15.10 29.43 低 度 
最 大 3 日 299 312.7 4.57 43.81 中 度 
最 大 7 日 243.9 239.3 -1.90 43.81 中 度 
最 大 30 日 153.7 171.9 11.82 28.09 低 度 
最 大 90 日 114.1 135.3 18.51 85.62 高 度 
零 流 量 天 数 0 9 76.92 高 度 
基 流 指数 0.20 0.07 -64.89 85.62 高 度 
年 极端 值 发 生 最 小 1 日 流量 发 生 时 间 4 36 800 22.51 低 度 
时 间 最 大 1 日 流量 发 生 时 间 205 212 3.83 31.09 低 度 
高 、 低 流量 发 生 年 低 流 量 发 生 次 数 7 8 14.29 14.25 低 度 
频率 及 持续 时 年 低 流量 平均 历时 3.5 5 42.86 26.5 低 度 
年 高 流量 发 生 次 数 6 7 16.67 25.31 低 度 
年 高 流量 平均 历 5 4.5 -10.00 65.38 中 度 
流量 变化 改变 流量 平均 增加 率 17 2.5 47.06 73.54 高 度 
率 及 频率 /次 流量 平均 减少 率 -2 -3.1 55.00 85.62 高 度 
每 年 流量 逆转 次 数 143 150 4.90 24.04 低 度 
流量 发 生 时 间 为 罗马 日 ,高 低 流量 发 生 频 率 及 持续 时 间 单 位 分 别 为 次 -年 "和 日 ,流量 变化 改变 率 单位 为 ms 


注 :流量 单位 为 ms ,年 极端 


,每 年 流量 逆转 次 数 为 次 .年 - 


(2) 年 极端 流量 。 变 化 后 ,最 小 1 日 .3 日 .7 日 
流量 均 呈 减 小 趋势 , 且 减 小 幅度 均 超 过 50% ;而 最 


大 流量 中 , 除 最 大 7 日 流量 外 , 均 呈 增加 趋势 ,其 中 
最 大 90 日 流量 增加 最 大 ,变化 率 为 18.51% ; 断 流 天 
数 从 变化 前 的 0 日 变 为 9 日 。 年 极端 流量 的 改变 度 


中 ,最 小 30 日 流量 .最 大 1 日 流量 和 最 大 30 日 流量 


分 别 为 28.09% .29.43% 和 28.09% ,属于 低 度 改变 ; 
最 小 90 日 流量 .最 大 3 日 流量 和 最 大 7 日 流量 的 改 


变 度 分 别 为 56.86% 、43.81% 和 43.81% ,为 中 度 改变 ; 
其 余 水 文 指标 的 改变 度 均 大 于 67% ,属于 高 度 改变 。 

(3) 年 极端 流量 发 生 时 间 。 变 化 后 ,最 小 1 日 
流量 和 最 大 1 日 流量 出 现 的 时 间 均 推迟 ,其 中 最 小 1 
日 流量 出 现 日 期 推迟 了 32 日 ;最 大 1 日 流量 的 发 生 
日 期 推迟 了 7 日 。 最 小 1 日 和 最 大 1 日 流量 的 改变 
度 分 别 为 25.51% 和 31.09%, 属 于 低 度 改 变 。 

(4) 高 . 低 流量 的 频率 及 历时 。 变 化 后 , 低 流量 
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图 4 黑河 上 游 营 落 峡 IHA 指标 改变 度 
Fig.4 The alteration degree of IHA index at Yingluoxia hydrologic station 


出 现 的 次 数 和 平均 历时 分 别 增长 了 14.29% 和 
42.86% ;高 流量 出 现 的 次 数 增加 了 16.67% ,但 高 流 
量 平 均 历时 减少 约 10%。 低 流量 发 生 次 数 .平均 历 
时 和 高 流量 发 生 次 数 的 改变 度 分 别 为 14.25% 、 
26.50% Fil 25.31% , 均 属 于 低 度 改变 ;而 高 流量 平均 
历时 的 改变 度 分 别 为 65.38% , 属于 中 度 改 变 。 

(5) 流量 变化 改变 率 及 频率 。 变 化 后 ,流量 平 
均 增 加 率 , 减 小 率 以 及 道 转 次 数 均 有 所 增加 ,其 中 
前 两 者 的 增加 幅度 在 50% 左 右 ,而 逆转 次 数 的 增加 
幅度 较 小 ,为 4.9%。 每 年 流量 逆转 次 数 的 改变 度 为 
24.04% , 属于 低 度 改变 ;而 流量 平均 增加 率 和 减 小 
率 发 生 显著 变化 ,属于 高 度 变 化 。 

从 表 1 中 可 以 看 出 ,黑河 上 游 年 内 径流 量 主要 
集中 在 6 一 10 月 ;1960 一 2015 年 黑河 上 游 月 平均 径 
流 呈 增加 趋势 ,年 极端 径流 中 多 个 指标 也 呈 增 加 趋 
热 。 王 杰 ™ 采 用 SWAT 模 型 预测 和 分 析 了 黑河 上 游 
径流 对 未 来 气候 变化 的 响应 ,结果 表明 ,年 内 径流 
主要 集中 在 6 一 9 月 ,上 且 研究 区 内 三 个 水 文 站 的 年 径 
流量 呈 增 加 趋势 ; 李 秋 菊 等 “采用 Budyko 假 设 下 的 
张 榴 经 验 公式 定量 分 析 了 黑河 上 游 径 流 的 变化 , 指 
出 1960 一 2015 年 营 落 峡 站 径流 量 显 著 增加 ,这 些 人 研 
究 验 证 了 本 文 的 分 析 。 

3.2 黑河 上 游 水 文 情势 改变 度 

进一步 分 析 黑 河上 游 径流 典型 年 份 的 变化 情 
况 ( 表 2), 可 以 发 现 ,1960 一 2002 年 黑河 上 游 以 平水 
年 和 枯 水 年 为 主 ,两 者 分 别 占 总 年 数 的 58.1% 和 
32.6% ;而 2003 一 2015 年 则 以 丰 水 年 为 主 , 占 总 年 数 


表 2 1960 一 2015 年 典型 年 份 的 分 布 情况 
Tab.2 Distribution of hydrological year types between 
1960 to 2015 


IN BE 统计 量 ”村 水 年 平水 年 丰 水 年 ED S 
1960 一 2002 年 “年数 14 25 4 
比例 /% 32.6 58.1 9.3 

2003 一 2015 年 ”年数 0 3 jo 216- 044 
比例 /% 0 23.1 76.9 


H 


注 :ED 表示 欧 拉 距 离 ;S 表 示 标 准 化 后 的 欧 氏 距离 ,数值 越 大 表示 


变化 越 大 。 下 同 。 


的 76.9% ,表明 黑河 上 游 的 径流 在 突变 年 份 前 后 两 
个 时 期 (1960 一 2002 年 与 2003 一 2015 年 ) 典 型 年 份 
组 成 发 生 了 显著 变化 。 同 时 分 析 突 变 年 份 前 后 各 
月 月 均 流 量 值 变化 情况 ( 表 1) ,可 以 看 出 年 内 6 个 月 
(3 月 .4 月 .8 月 ,9 月 10 月 11 月 ) 的 改变 程度 达到 
了 高 度 改 变 。 

根据 公式 (1) (2) 及 (3) 计 算 IHA 指标 的 改变 
FE ,采用 RVA 方 法 计算 的 整体 改变 度 为 54.68% 。 
根据 公式 (7) 和 (8) ,改进 的 RAV 方 法 计算 的 整体 改 
变 度 为 74.62%( 如 表 3 所 示 ) ,大 于 67% ,属于 高 度 改 
变 。 说明 RVA 方 法 计算 黑河 上 游 水 文 情势 变化 忽 
视 了 径流 典型 年 份 的 变化 ,造成 了 对 黑河 流域 水 文 
情势 整体 改变 度 的 低估 。 
3.3 黑河 上 游 水 文 情势 变化 原因 分 析 

从 上 文 的 分 析 结 果 可 以 看 出 ,黑河 上 游 河流 水 
文 情 势 发 生 了 高 度 改 变 。 水 文 情势 变化 的 原因 可 
分 为 人 类 活动 和 气候 变化 。 黑 河上 游人 类 活动 的 
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表 3 黑河 上 游 水 文 情势 整体 改变 度 计 算 结 果 
Tab.3 The alteration degree of the hydrological regime 


considering the characteristics of hydrological year types 


RVA 改进 的 RVA 
Dy 54.68% - 
S - 0.44 
OD - 74.62% 
改变 程度 中 高 


注 : 有 表示 RVA 法 计算 的 整体 改变 度 ;0D 表 示 由 改进 的 RVA 法 计 
算 的 整体 改变 度 。 


影响 微弱 3, 因此 ,本 文 对 黑河 上 游 周 边 4 个 气象 
站 的 观测 降水 和 气温 变化 进行 了 分 析 。 

通过 M-K 检 验 分 析 可 得 ,1960 一 2015 年 4 个 气 
象 站 的 观测 降水 均 呈 增加 趋势 ;4 个 气象 站 中 , FER 
站 和 野牛 沟 站 1960 一 2015 年 的 观测 降水 的 突变 年 
均 为 2002 年 ,而 张掖 站 和 祁 连 站 观测 降水 包含 多 个 


140, (a) 
120 2003 一 2015 年 
m1960 一 2002 年 
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图 5 2002 年 前 后 月 均 降 雨量 和 月 均 流 量 


突变 点 。 根 据 泰 森 多 边 形 法 加 权 计 算 的 黑河 上 游 
的 降水 量 ,2002 年 前 后 年 均 降 水 量 分 别 为 443.99 
mm 和 507.17 mm, 降水 增加 了 14.23%( 图 $)。 从 图 
5 可 以 看 出 ,降水 主要 集中 在 夏季 和 秋季 ,同时 , 夏 、 
秋季 也 是 降水 增加 量 较 大 的 季节 Ee 
研究 结论 一 致 ;除了 3 月 和 12 月 份 的 降水 呈 微 弱 减 
少 外 ,其 他 各 月 的 平均 流量 均 呈 不 同 程度 的 增加 ， 
增加 量 最 大 的 为 8 月 和 9 月 ,分别 增 加 14.8 mm 和 
21.96 mm。 降 水 的 增加 导致 平均 流量 改变 度 显著 变 
化 ,例如 8 月 .9 月 10 月 和 11 月 的 平均 流量 的 改变 
度 均 大 于 67% ,属于 高 度 改变 ,而 且 降 水 的 增加 还 
导致 其 他 极 值 流 量 成 不 同 程度 的 增加 ,如 最 大 90 日 
流量 的 改变 度 为 85.62% ,属于 高 度 改 变 (图 5)。 

4 个 气象 站 1960 一 2015 年 观测 气温 的 年 均值 如 
图 6 所 示 。4 个 气象 站 观测 温度 的 年 均值 均 呈 增加 
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Fig. 5 Monthly precipitation and runoff of upper Heihe River Basin from 1960 to 2015 


A 托 勒 站 温度 o 野牛 沟 站 温度 目 张 掖 站 温度 a 祁 连 站 温度 a "a z 
a = 
8 » ee 5 
Ba'n ee nM 0 0366x-65.285 
6 a a a R°=0.63, P=0.05 
4 
¥=0.0309x-60.316 
a R°=0.61, P=0.05 nt see Shr Se 
Y Aaa AA AAA Gyt rae Vie alat Ge aa 
o bea a aa a4 y=0.0374x-76.881 
4 R°=0.65, P=0.05 
A AwA-AAA 
2 BS na aang adSASe.S Ro o 
$040 40 0 a 5066 Goo” Oo he 
atas 6 6 aA y=0.0285x-59.521 
O 


R?=0.46, P=0.05 


1960 1970 1980 


2000 2010 2020 


图 6 黑河 上 游 周边 气象 站 年 均 气 温 及 变化 趋势 


Fig. 6 Mean annual temperature of meteorological stations in and around upper Heihe River basin 
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趋势 ,这 与 李 占 玲 等 ”的 研究 结果 一 致 。 一 方面 ， 
气温 的 升 高 增加 了 流域 蒸 散 发 ,造成 径流 的 减少 ， 
李 秋 菊 等 "研究 表明 ,2004 一 2015 年 , 莱 散 发 对 径 
流 增加 的 贡献 率 有 -35%; 另 一 方面 ,气温 的 升 高 导 
致 冰川 融 水 的 增多 ,导致 径流 的 增加 ,黑河 上 游 冰 
川 及 融雪 径流 约 占 总 径流 的 10%50 ,Sang 等 "研究 
结果 表明 温度 的 升 高 有 利于 径流 的 增加 。 王 宇 涵 
等 中 研究 表明 ,气温 的 升 高 对 径流 的 增加 有 贡献 ， 
这 是 因为 气温 升 高 导致 高 寒山 区 冻 土 活动 层 增 
厚 ,增加 了 土壤 蓄 水 容量 ,从 而 导致 降雨 下 渗 量 增 
加 和 基 流 量 增 大 。 总 体 而 言 ,黑河 上 游 对 气候 变 
化 敏感 ', 黑 河上 游 径流 的 增加 主要 是 受 “ 暧 湿 " 气 
候 影 响 呈 ,降水 与 气温 升 高 均 有 助 于 径流 的 增加 ， 
导致 春季 径流 增加 的 主要 因素 为 温度 ,降水 次 之 ; 
而 导致 夏秋 、 秋 季 径 流 增 加 的 主要 因素 为 降水 , 温 
BR, 


14S 


4 结论 


=A 

径流 典型 年 份 ( 丰 、 平 .村 水 年 ) 的 变化 对 河流 
生态 系统 的 结构 和 功能 具有 重要 影响 。 改 进 的 
RVA 法 将 径流 典型 年 份 的 变化 纳入 整体 改变 度 的 
计算 过 程 中 ,相对 全 面 地 评价 径流 水 文 情势 的 变 
化 。 改 进 RVA 与 原来 的 RVA 对 数据 的 要 求 相似 ， 
具有 普遍 适用 性 。 本 文采 用 改进 的 RVA 分 析 了 黑 
河上 游 河 流水 文 情势 的 变化 情况 。 整 体 而 言 ,黑河 
上 游 河 流水 文 情势 发 生 了 高 度 改变 ,黑河 上 游 各 月 
平均 径流 年 极端 流量 大 多 数 指标 和 流量 变化 改变 
率 以 及 逆转 次 数 基 本 呈 不 同 程度 的 增加 。 黑 河上 
游 水 文 情 势 的 变化 与 降水 等 气候 因素 密切 相关 , 气 
修 变 化 是 导致 黑河 上 游 水 文 情势 高 度 改变 的 主要 
原因 。 
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Flow regime alterations of upper Heihe River based on improved RVA 
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and Engineering, Southern University of Science and Technology, Shenzhen 518055, Guangdong, China; 


3. College of River and Ocean Engineering, Chongqing Jiaotong University, Chongqing 400074, China) 


Abstract: Objective evaluation of flow regime alteration is the foundation of river ecological restoration and 
integrated water resources management. Various approaches that take hydrological year types into consideration 
could objectively evaluate the altered flow regime. Due to climate change in the recent decades, the flow regime 
of upper Heihe River has altered and needs further research. Our study area was the upper Heihe River, where we 
employed an improved range of variability approaches to analyze the alterations of the flow regime. Firstly, we 
applied the Mann-Kendall method to identify the catastrophe point. We analyzed the flow regime using 
hydrological alteration indicator method. This method includes 33 parameters classified into five groups: The 
magnitude of monthly water conditions, magnitude and duration of annual extreme water conditions, timing of 
annual extreme water conditions, frequency and duration of high and low pulses, and rate and frequency of water 
condition changes. Based on these groups, we used the range variability approach to calculate the degree of 
alteration of these 33 parameters. We measured the difference in typical hydrologic years before and after the 
catastrophe point using the Euclidean distance. The degree of hydrologic alteration was calculated considering 
the alteration of these 33 parameters and the difference between typical hydrologic years. Our results indicated 
that the year 2002 was catastrophic point of the river flow regime. The main hydrological type years before and 
after 2002 were normal flow and wet years. The improved method calculated an alteration degree of 74.62%, 
while the range of variability approach showed an alteration of 54.68%, indicating underestimation of the effect. 
Compared to 1960-2002, the monthly mean streamflow in 2003-2015 increased to various extents in the upper 
Heihe River’ s flow regime. The annual minima 1 day, 3 day, 7 day, and annual maxima 7 day decreased, whereas 
the other extreme water condition indicators increased. The frequency and duration of high and low pulses 
changed slightly, and the average peak flow duration decreased while the other indicators increased. Rate and 
frequency of water condition increased. The degree of hydrologic alteration was 74.62% considering the change 
of HYTs, which indicated a high alteration in the upper Heihe River’ s flow regime. Climate change is the main 
reason causing the flow regime alteration. 

Keywords: alteration of river flow regime; range of variability approach; hydrological year types; upper Heihe 


River 


